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We study how the neutron-star equation of state affects 
the onset of the dynamical instabihty in the equations of mo- 
tion for inspirahng neutron-star binaries near coalescence. A 
combination of relativistic effects and Newtonian tidal effects 
cause the stars to begin their final, rapid, and dynamically- 
unstable plunge to merger when the stars are still well sepa- 
rated and the orbital frequency is « 500 cycles/sec {i.e. the 
gravitational wave frequency is approximately 1000 Hz). The 
orbital frequency at which the dynamical instability occurs 
(i.e. the orbital frequency at the innermost stable circular or- 
bit) shows modest sensitivity to the neutron-star equation of 
state (particularly the mass-radius ratio, M/Ro, of the stars). 
This suggests that information about the equation of state of 
nuclear matter is encoded in the gravitational waves emitted 
just prior to the merger. 



PACS Numbers: 97.80. Fk, 04.25.Dm, 04.40. Dg, 97.60. Jd 



I. INTRODUCTION 

Binary neutron star systems which are spiraling toward 
their final coalescence under the dissipative influence of 
gravitational radiation reaction forces are the primary 
targets for detection of gravitational waves by interfer- 
ometric gravitational wave detectors such as LIGO and 
VIRGO Extracting the gravitational waves from 

the detector noise and making use of the information 
encoded in the signals will require a thorough knowl- 
edge of the expected waveforms produced by these bi- 
naries In this paper we explore the effect of the 
neutron-star equation of state on the orbital evolution 
and gravitational wave emission of binaries just prior to 
merging. Specifically, we show that a combination of 
post-Newtonian (relativistic) effects and Newtonian tidal 
effects (which depend on the equation of state) conspire 
to induce a dynamical instability in the orbital motion, 
which causes the plunge to final coalescence to begin 
somewhat sooner - and to proceed somewhat faster - 
than it would simply under the influence of the dissi- 
pative radiation reaction force. Thus the motion of the 
bodies during the late stages of binary inspiral depends 
on the structure of the neutron stars. Consequently, the 



gravitational waveform emitted during this short portion 
of the final coalescence will be imprinted with informa- 
tion about the nuclear equation of state. 

During the final ^10 minutes or the last ~ 8000 orbits 
of a neutron star binary inspiral, the orbital frequency 
increases from about 5 Hz on up to a cutoff of a few 
hundreds to a thousand Hertz (roughly corresponding 
to the orbital frequency when the final plunge begins). 
Thus the gravitational wave frequency (twice the orbital 
frequency for the dominant quadrupole radiation) chirps 
through the LIGO detector bandwidth during this period 

. The evolution of the binary in these last few minutes 
of the inspiral is very sensitive to a number of relativis- 
tic effects, such as gravitational-wave tails and spin-orbit 
coupling (dragging of inertial frames) . The gravitational 
waveform emitted by the binary in this portion of the 
coalescence, the adiabatic inspiral, is currently being ex- 
tensively studied . During most of this inspiral phase 
the neutron stars can be treated as simple point masses 
because the effects associated with the finite stellar size 
turn out to be small: (i) The neutron star has too small 
a viscosity to allow for angular momentum transfer from 
the orbit to the stellar spin via viscous tidal torque ; 
(ii) The effect of the spin-induced quadrupole is negli- 
gible unless the neutron star has rotation rate close to 
the break-up limit (ui) Resonant excitations of 

neutron star internal modes (which occur at orbital fre- 
quencies less than 100 Hz) produce only a small change 
in the orbital phase due to the weak coupling between 
the modes and the tidal potential pdj-p^; (iv) The cor- 
rection to the equation of motion from the (static) tidal 
interaction is of order {Ro/r f' (where Ro is the neutron 
star radius, r is the orbital separation jlj]), which is neg- 
ligible except when r is smaller than a few stellar radii. 
Since Ro ~ 5Af for a typical neutron star of mass 1.4Mq 
and radius 10 km, the tidal effect is essentially a (post)^- 
Newtonian correction [||. The expression for the phase 
error induced by the tidal effect is given in Ref. ||l^ . The 
fact that the evolution of the binary system as it sweeps 
through the low frequency band of the detector is insen- 
sitive to finite-size effects means that the measurement 
of the inspiral waveform will allow us to probe cleanly 
into the intricate structure of general relativity, and to 
test whether general relativity is the correct theory of 
gravity |l^Jl^. Moreover, some of the parameters of the 
binary system, such as the masses of the stars, can be 
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determined with reasonable accuracy [ |l7[|lq ]. However, 
the waveform's lack of dependence on the finite size of 
the objects during the most of the adiabatic inspiral also 
implies that information about the internal structure of 
the neutron star is only imprinted on the radiation emit- 
ted just prior to coalescence when the orbital radius is 
small. 

Indeed at small orbital separations, tidal effects are 
expected to be very important. In a purely Newtonian 
analysis, the interaction potential between star M' and 
the tide-induced quadrupole of M, Vude ~ —M''^R^^/r^, 
increases with decreasing r. The potential becomes so 
steep that a dynamical instability develops, accelerating 
the coalescence at small orbital radius. This Newtonian 
instability has been fully explored using semi-analytic 
models in Ref. and Ref. (hereafter referred to 
as LS). It has also been examined numerically in Refs. 

However, a purely Newtonian treatment of the binary 
at small separation is clearly not adequate, as general 
relativistic effects will also be important in this regime; 
and general relativistic effects can also make the orbit 
unstable. For example, a test particle in circular orbit 
around a Schwarzschild black hole will experience an "in- 
nermost stable circular orbit" at risco = 6Af (or hM in 
harmonic coordinates). This unstable behavior is caused 
by higher-order relativistic corrections included in the 
Schwarzschild geodesic equations of motion. For comput- 
ing the orbital evolution of two neutron stars of compara- 
ble mass near coalescence, the test-mass limit is obviously 
inadequate. In order to explore the orbital instability 
for such systems, Kidder, Will and Wiseman ||2^ (here- 
after referred as KWW) developed hybrid equations of 
motion. These equations augment the the Schwarzschild 
geodesic equations of motion with the finite-mass terms 
of the (post)^/^-Newtonian equations of motion. Includ- 
ing these finite mass terms in the equation of motion 
moved the innermost-stable-circular-orbit radius farther 
out (in units of the total mass).|^ In this paper, we aug- 
ment the hybrid equations with contributions due to the 
tidal deformation of the stars. In a nutshell, the work 
presented here combines the Newtonian tidal analysis of 
LS with the relativistic point-mass analysis of KWW 
p2| to yield a more complete picture of the neutron-star 
coalescence prior to merging. 

We note, that unlike a test-particle around a 
Schwarzschild black hole, the very notion of "innermost 
stable circular orbit" is poorly defined for objects of com- 
parable mass. After all, in the relativistic regime the bi- 
nary orbit will be decaying rapidly due to radiation reac- 



^See Wex and Schafer for a critique and an alterna- 
tive construction. Their post-Newtonian calculation suggests 
that the innermost stable orbit may occur at an even greater 
separation. 



tion; thus the orbit is not circular, but rather a decaying 
spiral. In order to give a semi-quantitative definition of 
"innermost stable circular orbit" we use the artifice of 
"shutting off" all the dissipative terms in the equation of 
motion and looking for the point where the solutions of 
the remaining non-dissipative equations become dynam- 
ically unstable. The use of hybrid equations of motion 
augmented with the tidal terms allows us to map out the 
dependence of the critical radius Hsco, or the correspond- 
ing orbital frequency /isco, for a wide range of allowed 
neutron star equations of state (parametrized by radius 
and effective polytropic index; see Figure 1). We believe 
that clearing up such dependence is important, and this 
analysis provides a benchmark with which comparisons 
can be made with future numerical results. Indeed, an 
important point we wish to make in this paper is that nei- 
ther relativistic (post Newtonian) effects nor Newtonian 
tidal effects can he neglected near the instability limit, and 
the critical frequency can be much lower than the value 
obtained when only one of these effects are included. 
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FIG. 1. The critical orbital frequency (at the inner-most 
stable orbit) as a function of the ratio of the neutron star 
radius Ro and mass M. The solid curves show the results 
including both relativistic and tidal effects (the lower curve is 
for r = 3 while the upper one is for F = 2) , the dashed curves 
are the Newtonian limit given by Eqs. (8)-(9). The vertical 
line corresponds to Ro/M = 9/4, the minimum value for any 
physical neutron star. The insert is a close-up for the nominal 
range of Ro/M = 4 — 8 as given by all the available nuclear 
EOS's. Two curves within the insert should bracket all the 
physical values of /isco. 

The main results of our analysis are summarized in 
Figure 2 and 3. Figure 2 shows that the rate of radial 
infall for stars near coalescence is substantially underes- 
timated if one models the coalescence as a Newtonian 
circular-orbit decaying solely under the influence of ra- 
diation reaction (top dotted curve). In other words, the 
rate of coordinate infall is substantially enhanced by the 
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non-dissipative terms. Somewhat more relevant for ob- 
servational purposes, Figure 3 shows that the number of 
orbits (or gravitational wave cycles) per logarithmic fre- 
quency interval is substantially reduced by the unstable 
collapse of the orbit. Both plots show modest sensitivity 
to the equation of state. 
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FIG. 2. The radial infall coordinate velocity during bi- 
nary coalescence, with M = M' = 1.4Mq, Ro/M = 5, F = 3, 
all calculated using 2.5PN radiation reaction. The solid 
line is the result including relativistic and tidal effects, the 
short-dashed line includes only tidal effects, the long-dashed 
line includes only relativistic effects. The dotted line is the 
point mass "Newtonian" result. 
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FIG. 3. The number of orbits the binary spends per log- 
arithmic frequency. The labels are the same as in Fig. 2. 



The remainder of the paper is organized as follows: In 
section II we present our equations of motion. In sec- 
tion III we examine the orbital stability using the non- 
dissipative portion of the equations of motion, and thus 
identify the location of the "innermost stable circular or- 
bit" . In section IV we include the dissipative terms that 
were omitted in the analysis of section III, and evolve the 
full equations of motion. In section V we briefly discuss 
the relevance of our results to numerical hydrodynamic 
calculations and to gravitational wave signal analysis. 



II. EQUATIONS OF MOTION INCLUDING 
TIDAL AND GENERAL RELATIVISTIC 
EFFECTS 

Consider a binary containing two neutron stars of mass 
M and M' , each obeying a polytropic equation of state 
P = Kp^. We use the compressible ellipsoid model for 
binary stars developed in LS . Basically, we model the 
tidally deformed neutron star as an ellipsoid, with inter- 
nal density profile similar to that of a spherical polytrope. 
The dynamics of such a neutron star (so called Riemann- 
S ellipsoid) is characterized by the three principal axes 
(fli, a2, as for star M and 0^,02,03 for star M'), the an- 
gular velocity (Jl and of the ellipsoidal figure about 
a principal axis (perpendicular to the orbital plane) and 
the internal motion of the fluid with uniform vorticity. 
The non-zero internal fluid motion is necessary because 
the binary neutron stars are not expected to corotate 
with the orbit due to rapid orbital decay and small vis- 
cosity 0,^. Although the Newtonian tidal interaction 
between the neutron stars can be treated exactly in the 
linear regime using mode decomposition the ellip- 
soid model has the advantage that it can be extended 
to the nonlinear regime at small orbital radii, when the 
tidal deformation of the star becomes significant. 

The Newtonian dynamical equations for the binary 
neutron stars as derived in LS include the familiar New- 
tonian (l/r^) force-law for point-masses orbiting one an- 
other; they further contain Newtonian terms involving 
finite size (tidal) effects. A post-Newtonian treatment 
of the tidal problem would give the relativistic correc- 
tions to these terms, namely the standard point-mass, 
post-Newtonian corrections to the equations of motion, 
as well as relativistic corrections to the quadrupole mo- 
ment and corrections due to higher moments the bodies. 
(See Appendix F of ||2^.) To insure that our equations 
of motion at least agree with the known post-Newtonian, 
point-mass equations we augment these Newtonian equa- 
tions of motion with the hybrid equations of KWW. How- 
ever, we use only the Newtonian equations to describe 
the evolution of the neutron stars structure (a^ and 
and the fiuid motion (the figure rotation rate and the 
internal vorticity) within the stars. These are given by 
Eqs. (2.18)-(2.22) of LS. In other words, we neglect the 
relativistic corrections to the fluid motion, self-gravity 
and tidal interaction. These corrections are secondary 



3 



effects and should not modify the orbital dynamics ap- 
preciably (e.g., the Newtonian tidal interaction between 
the two stars scales approximately as M'^R^/r^, and its 
relativistic correction is of order M/r smaller). As noted 
before, the tidal interaction enters the Newtonian po- 
tential as a correction of 0[{Ro/r)^] ^ 0[{AIt/r)^], in 
effect, as a (post)^-Newtonian term. Therefore, by not 
including relativistic corrections to these tidal terms, we 
are merely omiting terms which are of (post)^-Newtonian 
order. In fact, the largest error comes from neglecting the 
post-Newtonian correction (of order M/Ro) to the inter- 
nal stellar structure (see Sec. III.B for an estimate of its 
effect on /isco)- The relativistic corrections to the orbital 
motion, however, are very important. Our equations of 
orbital motion can be assembled from Eqs. (2.23)-(2.24) 
of LS and from Eqs. (1.2)-(1.3) of KWW: 
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where Mt = M + M' is the total mass, a (a') is the mis- 
alignment angle between the tidal bulge of M {M') and 
the line joing the two dimensionless 
structure constants depending on the mass concentra- 
tion within the stars. In Eqs. (1) and (2) the last two 
lines contain the "dissipative" terms due to gravitational 
radiation reaction. The quantities Ah, Bh, ^5/21 ^5/21 
which include the "hybrid" corrections to the equation of 
motion, are given by 
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where = P + r^0^, 77 = fi/Mt and /i MM' /Mt. Also 
the multipoles moments can be expressed as 



K„Ma2/5, 4 = <M'af /5 



(7) 



Note that in Eqs. (l)-(2), we have also included the lead- 
ing order radiation reaction forces due to tidal deforma- 
tion. 

Admittedly, this is not a consistent post-Newtonian 
expansion of the true equations of motion; however it 
is correct in several important limiting cases: (i) In the 
limit that a^ — > and — > and the limit 77 ^ 0, 
we recover the exact Schwarzschild equation of motion, 
(ii) In the point-mass limit (a^ —t and a[ 0) we 
recover the hybrid equations given in KWW. KWW pre- 
sented an argument that suggested that the higher-order, 
r7-dependent, (post)'^-Newtonian - as yet uncalculated - 
corrections to these equations have only a modest effect 
on the equations of motion. See Figure 6 of Ref. [ p2[ . 
However, until these terms are calculated it is unclear 
just how large an effect they will have on the location 
of the innermost stable orbit, (iii) In the non- relativistic 
limit we recover the equations of motion given in LS. 
These equations contain the dominant contributions to 
the equations of motion due to the finite sizes of the ob- 
jects. 

Note that although Eqs. (1) and (2) make reference 
to the orbital radius r, we are always aware that this is 
a gauge dependent quantity and of little meaning for a 
distant observer. Observationally, the more meaningful 
quantity is the orbital frequency as measured by distant 
observers, and we shall use frequency rather the radius 
in presenting most of our results. 



III. INSTABILITY OF THE NON-DISIPATIVE 
EQUATIONS OF MOTION 

A. Method to Determine the Stability Limit 

We now form a set of non-dissipative equations of mo- 
tion, by simply discarding the gravitational radiation re- 
action terms given in the last two lines of Eq. (1) and 
Eq. (2). These non-dissipative dynamical equations ad- 
mit equilibrium solutions, which are obtained by setting 
r ~ = 9 — Qorb = a = a' = as well as hi = a[ = 0. 
For a given r, the evolution equations for the neutron 
star structure reduce to a set of algebraic equations for 
at and a[, while the orbital equation (2) gives the or- 
bital frequency ^orb- These equations are solved using a 
Newton-Raphson method, yielding an equilibrium binary 
model. Thus a sequence of binary models parametrized 
by r can be constructed. 
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To determine the stability of the orbit of a binary 
model, we simply use the equilibrium parameters as ini- 
tial conditions for our non-dissipative equations of mo- 
tion. We add a small perturbation to the equilibrium 
model and let the system evolve. In this way we locate 
the critical point of the dynamical equations, correspond- 
ing to the dynamical stability limit of the equilibrium bi- 
nary or the inner- most stable circular orbit: for r > risco, 
the binary is stable, and the system oscillates with small 
amplitude about the initial configuration; for r < risco, 
the binary is unstable, and the perturbation grows, lead- 
ing to the swift merger of the neutron stars even in the 
absence of dissipation. 

B. Results 

For concreteness, we present results only for binary 
neutron stars with equal masses {M = M'), both hav- 
ing zero spin at large orbital separation, although our 
equations are adequate to treat the most general cases 

The polytropic relation P — Kp provides a use- 
ful parametrization to the most general realistic nuclear 
equation of state (EOS). Since the radius Ro of the non- 
rotating neutron star of mass M is uniquely determined 
by K and F, we can alternatively use Ro/M and F to 
characterize the EOS. For a canonical neutron star with 
mass M = 1.4Mq, aU EOS tabulated in ^ give Ro/M 
in the range of 4 — 8, while modern microscopic nuclear 
calculations typically give Ro/M = 5 For a given 

Ro/M ^ the polytropic index F specifies the mass concen- 
tration within the star. Except for extreme neutron star 
masses (M <, O.SM© or M ^ I.8M0) typical values of F 
lie in the range of F = 2 - 3 

In Table I, we list the physical properties of the equi- 
librium binary neutron stars at the dynamical stability 
limit for several values of Ro/M and F = 3. In Figure 1, 
the orbital frequency /isco is shown as a function of Ro/M 
for F = 2 and F = 3. Clearly, in the limit of Ro/M 0, 
/isco approachs the point mass result /isco = 697M^l 
Hz obtained in KWW0. In the non-relativistic limit wc 
recover the pure Newtonian result [ p^|jis[| : 

/isco = 657Mi:l(5M/i?o)^/2 (jjz) (r = 3), (8) 
/isco = imM^libM / Rof'^ (Hz) (F = 2). (9) 

For typical neutron star radius Ro/M — 5, the critical 
frequency ranges from 488 Hz (for T — Z) to 540 Hz (for 
T = 2), while both the pure Newtonian (with tides) cal- 
culation and the pure point-mass hybrid calculation give 
a result 30 — 40% larger. There are two physical causes 



^KWW give a correct expression for rism/Mt, but in- 
correctly give /isco=710 Hz due to a numerical error. 



for the reduction in /isco: (i) The binary becomes unsta- 
ble at larger orbital separation due to the steepening of 
the interaction potential from both tidal and relativistic 
effects; (ii) For a given orbital radius (itself a gauge de- 
pendent quantity) , the post-Newtonian orbital frequency 
as measured by an observer at infinity is smaller than the 
Newtonian orbital frequency^. We conclude that to ne- 
glect either the tidal effects or the relativistic effects can 
lead to large error in the estimated critical frequency. 

Except for the intrinsic uncertainties associated with 
the hybrid equations of motion p^ , p3| , the main uncer- 
tainty in our determination of /isco comes from neglect- 
ing post-Newtonian corrections to (i) the stellar struc- 
ture and (ii) the the tidal potential. The first correc- 
tion decreases the tide- induced quadrupole; the fractional 
change is of order —M/Ro- The second increases the 
quadrupole by a fraction of order M' /r. We can estimate 
how much /isco is modified by these two corrections. In 
Newtonian theory, risco is approximately determined by 
the condition MM' /r ~ M''^R\/r^. Including the rel- 
ativistic corrections this condition becomes MM' /r ~ 
Af'2i?5/r6(i - 5)^ where 5 ~ [0{M / Ro) - 0{M' /r)] <, 
20%. Thus the change in /isco due to these two effects 
is A/isco//isco — 0.3(5 < 6%, i.e., the critical frequency 
increases by a few percent |^ . 

As emphasized in Sec. I, the critical radius (or crit- 
ical frequency), at which the non-dissipative equations 
become dynamically unstable, is meaningful only in the 
sense that when r < risco, the binary will coalesce on 
dynamical (orbital) timescale even in the absence of dis- 
sipation. In the realistic situation, the dissipative radia- 
tion reaction forces will also be rapidly driving the binary 
to coalescence. Therefore to determine the significance 
of the dynamical instability we must compute the orbital 
evolution with the full equations of motion - including 
the radiation reaction. 



IV. ORBITAL EVOLUTION PRIOR TO MERGER 

We now include the dissipative radiation reaction 
forces in our analysis. In this case the plunge will 
be driven by both the dissipative, as well as the non- 
dissipative effects associated with the steepening poten- 
tial (both the tidal potential and the relativistic poten- 
tial). But what effect is dominant? In order to numer- 
ically investigate this question, we choose a specific sys- 
tem with M = M' = IAMq, Ro/M = 5 and F = 3. The 
orbital evolution begins when the stars are well outside 
the innermost stable circular orbit limit. We consider 
four different inspiral scenarios, (i) A purely dissipa- 
tive inspiral: a system of point masses subject only to 



In the case of equal masses, at first post-Newtonian order 
^orb = ^Kepier[l " {11 / 8) (Mt / r)]. See Ref. §. 
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a Newtonian (l/r^) force and (post)^/^-Newtonian radi- 
ation reaction force. In this case the infaU rate is given by 
Vr = r = — (64/5)?7(Aff/r)^. [Specifically, we set Ah — 1 
and Bh ^ a, = a'^ = hk = I'uk = ^ in Eqs. (l)-(2).] 
This is depicted by the dotted curve in Figure 2 and 3. 
(ii) A purely relativistic plunge in which we neglect the 
tidal effects [Specifically we set = a- = Ikk = I'ku — 
in Eqs. (I)-(2)]. This relativistic case is depicted by the 
long-dashed curve in Figure 2 and 3. (iii) A tidally en- 
hanced plunge: we include only the Newtonian terms in 
Eqs. (f )-(2) and the radiation reaction force. [Specifically 
we set Ah — 1 and Bh = in Eqs. (I)-(2).] This case 
is depicted by the short-dashed curve in Figure 2 and 3. 
(iv) Finally, we evolve the complete dynamical equations 
including all terms in Eqs. (l)-(2); this is depicted by 
the solid curve in Figure 2 and 3. Each intergration is 
terminated when the surfaces of the stars touch, i.e., at 
r ~ 2ai (for the point-mass problem, the calculation is 
terminated at r ~ 2i?o). 

In Figure 2 we clearly see that the non-dissipative ef- 
fects - tidal and relativistic - substantially increase the 
rate of infall. The radial velocity at binary contact is 
comparable to the tangential velocity. We also note that 
the radial coordinate velocity is a gauge dependent quan- 
tity; therefore our only intent in using it in Figure 2 is 
to convey the general trend that the rate of infall is en- 
hanced by the dynamical instability. 

Figure 3 shows the number of orbits the binary spends 
per logarithmic frequency. In the simplest point-mass, 
Newtonian-plus-radiation- reaction case [case (i) above], 
the result can be calculated analytically 

dNorb/dlnforb = (5/1927r)^-iAf,^'/'(2^/o,f,)-5/3 

= 1.95 X 105(Afi.4/ort/Hz)-5/3^ (10) 

which gives 6 cycles at /^^ ~ 2/or& = 1000 Hz. In con- 
trast, the tidal and relativistic effects reduce this number 
to less than 2. 

Figure 4 shows the wave energy emitted around a given 
frequency, dEg^/dln forb = i^orb/^orb)Egu,, where Eg^ 
is calculated using the simple quadrupole radiation for- 
mula. The Newtonian plus radiation reaction result is 
dEg^/dlnforb = 1.63 X 10-3(/„rt/Hz)2/3MVi?„. We 
see that the radiation power near contact becomes much 
smaller. Note that dEgm/dlnforb calcuated in this way 
is not exactly the energy power spectrum, which must 
be obtained from the Fourier transform of the waveform 
[ p8[ ; however, it does provide a semi-quantitative fea- 
ture of the full analysis; in particular, the dip in the 
dEgw I d In /orb curve around 600 Hz results from the dy- 
namical instability of the orbit (see also Refs. [^^, al- 
though the calculations presented there are purely New- 
tonian) . 
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FIG. 4. The quadrupolar gravitational energy emitted 
near a given frequency. The labels are the same as in Fig. 2. 



V. DISCUSSION 

A number of authors have tried to define and locate the 
innermost stable circular orbit for relativistic coalescing 
systems of comparable masses p^j30| "|3^ in order to char- 
acterize the final moments of a binary coalescence (See 
for a discussion). The results of the various analy- 
ses are not converging to an agreed answer. Obviously, 
the precise nature of the final coalescence of two neutron 
stars will only be determined by a numerical simulation 
using full general-relativistic hydrodynamics. However, 
the present analysis does point to two intersting features 
to look for in a full numerical treatment: (i) To get even 
a qualitative picture of the coalescence, it is necessary to 
begin the numerical evolution when the stars are still well 
separated, i.e. before the onset of the orbital dynamical 
instability. The instability - the plunge - causes the coa- 
lescence to proceed much more swiftly than a coalescence 
driven solely by radiation reaction; thus the actual coa- 
lescence may differ qualitatively from one computed with 
a simple radiation-reaction driven inspiral. The final co- 
alescence may be more of a splat, than the slow winding 
together of the stars [^^ . (ii) The instability results from 
both the tidal effects and the relativistic corrections in 
the equations of motion. KWW showed that there is no 
instability in the first post-Newtonian relativistic equa- 
tions of motion; the instability does not show up until at 
least second post-Newtonian order. Therefore, in order 
for numerical simulations to see the effects of the rela- 
tivistic unstable plunge, it will probably require the use 
of at least second-order, post-Newtonian hydrodynamic 
code. KWW also showed that the location of the dy- 
namic instability does not converge very rapidly as one 
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increases the post- Newtonian order of the approximation. 
[This fact led KWW to the introduce the hybrid equa- 
tions of motion.] Thus, to get even a quahtatively ac- 
curate evolution of the binary near coalescence, it may 
be necessary to use a full general relativistic hydrody- 
namic treatment of the coalescence problem, and begin 
the evolution when the stars are still well separated. 

As we have shown, the dynamical instability in the 
equation of motion will, in effect, cut off the chirp- 
ing waveform. The frequency of the cut-off is some- 
what dependent on the neutron star equation of state. 
Only in this late stage of the evolution does the equa- 
tion of state leave a tell-tale sign in the emitted wave- 
form. However, devising a strategy to dig this infor- 
mation from the detector output requires further con- 
sideration. Most detection/measurement strategies for 
coalescing binaries involve integrating template wave- 
forms against long stretches (8000 orbits!) of raw out- 
put data, the idea being that one can detect/measure a 
relatively low amplitude signal by integrating for a long 
time. Looking for the signature of this very late stage 
plunge is precisely the opposite: we are looking at the 
waveform just before coalescence when the amplitude is 
fairly strong, but the plunge is of fairly short in dura- 
tion. So answering questions about the plunge (such as, 
at what orbital frequency did it begin?) requires mea- 
suring a relatively large amplitude, but short-duration, 
effect. Clearly, analysis of such events will require a dif- 
ferent detection strategy [ p8| . 
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TABLE L Physical quantities at the inner-most stable or- 
bit (dynamical stability limit) of neutron star binary, with 
M = M', r = 3, and zero spin at large orbital radii. Here Ro 
is the neutron star radius, ai, 02, 03 are the axes of the ellip- 
soidal neutron star (ai along the binary axis, 03 perpendicular 
to the orbital plane), and Mi. 4 = M/(1.4M0). 
The case in the top row is the purely Newtonian cal- 
culation using LS equations of motion. The last row is 
the point-mass calculation using the hybrid equations of 
motion of KWW. 
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FIG. 5. The critical orbital frequency (at the inner-most 
stable orbit) as a function of the ratio of the neutron star 
radius Ro and mass M. The solid curves show the results 
including both relativistic and tidal effects (the lower curve is 
for r = 3 while the upper one is for F = 2), the dashed curves 
are the Newtonian limit given by Eqs. (8)-(9). The vertical 
line corresponds to Ro/M — 9/4, the minimum value for any 
physical neutron star. The insert is a close-up for the nominal 
range of Ro/M = 4 — 8 as given by all the available nuclear 
EOS's. Two curves within the insert should bracket all the 
physical values of /iaco- 



FIG. 6. The radial infall coordinate velocity during binary 
coalescence, with M = M' = I.4M0, Ro/M = 5, F = 3, 
all calculated using 2.5PN radiation reaction. The solid 
line is the result including relativistic and tidal effects, the 
short-dashed line includes only tidal effects, the long-dashed 
line includes only relativistic effects. The dotted line is the 
point mass "Newtonian" result. 



FIG. 7. The number of orbits the binary spends per loga- 
rithmic frequency. The labels are the same as in Fig. 2. 



FIG. 8. The quadrupolar gravitational energy emitted near 
a given frequency. The labels are the same as in Fig. 2. 
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